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Abstract 

We  present  existence,  uniqueness  and  continuous  dependence  (with  respect  to  probability  distribu¬ 
tions  on  polarization  parameters)  of  solutions  in  Maxwell  systems.  This  provides  a  theoretical  and 
computational  foundation  for  associated  inverse  problems. 
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1  Introduction 

In  this  note  we  consider  well-posedness  questions  for  the  variational  solutions  of  one  dimensional  Maxwell’s 
equations  with  an  absorbing  left  boundary  condition,  a  supraconducting  right  boundary  condition  and  a 
general  macroscopic  polarization  term  which  includes  uncertainty  in  the  dielectric  parameters.  For  these 
solutions,  we  establish  existence,  uniqueness  and  continuous  dependence  on  the  uncertainty  measures  in 
a  Prohorov  metric  (see  [B],  [BB]  for  definitions  and  details)  sense.  As  explained  below,  these  results  can 
be  readily  used  in  an  inverse  problem  methodology  to  determine  the  unknown  distribution  of  the  dielec¬ 
tric  parameters  which  govern  the  behavior  of  the  electric  held  and  the  electric  polarization  in  a  general 
heterogeneous  material  with  multiple  mechanisms  (Debye,  Lorentz,  etc.)  and  relaxation  parameters. 


2  Problem  Formulation 


We  consider  the  one  dimensional  problem  formulation  in  [BBL].  Assuming  D  =  eE  +  P,  we  may  obtain 
Maxwell’s  equations  [E]  in  second  order  form  given  by: 

HoeE  +  fi,0InP  +  Ho —  E"  =  —  Ho Js  in  U  fl0,  (1) 


where  E  is  the  transverse  component  of  the  electric  held,  P  is  the  material  macroscopic  electric  polar¬ 
ization,  e  =  e(z)  is  the  dielectric  permittivity  and  a  =  <j(z)  is  the  conductivity  of  the  material.  The 
boundary  conditions  that  we  are  assuming  are  absorbing  at  z  =  0  and  supraconducting  at  z  =  1: 


E  -  cE'  1  =  0,  E(t,  1)  =  0. 

J  2  =  0 


(2) 


Our  initial  conditions  are 

E(0,z)  =  *(z),  E{Q,z)  =  ^(z).  (3) 

To  describe  the  behavior  of  the  electric  polarization  P,  we  begin  with  the  general  formulation  of 
Chapter  2  of  [BBL]  by  employing  a  polarization  kernel  g  in  the  convolution  expression 

rt 

(4) 


P(t,z)=  /  g(t  —  s,  z’,t)E(s,  z)ds. 

_  Jo 
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As  explained  in  [BBL],  this  general  formulation  includes  as  special  cases  the  well  known  orientational  or 
Debye  polarization  model,  the  electronic  or  Lorentz  polarization  model,  and  linear  combinations  thereof, 
as  well  as  other  higher  order  models.  In  the  Debye  case  the  kernel  is  given  by 

g(t;  r)  =  e0(es  -  £oo )/r  e~t/r, 

while  in  the  Lorentz  model  (again  see  [BBL]),  it  takes  the  form 

g(t;r)  =  toLOp/vo  e~t/2r  sin(t'ot). 


However,  use  of  these  kernels  presupposes  that  the  material  may  be  sufficiently  defined  by  a  single 
relaxation  parameter  r,  which  is  generally  not  the  case.  In  order  to  account  for  multiple  relaxation 
parameters  in  the  polarization  mechanisms,  we  allow  for  a  distribution  of  relaxation  parameters  which 
is  conveniently  described  in  terms  of  a  probability  measure  F.  Thus,  we  define  our  polarization  model 
in  terms  of  a  convolution  operator 


Q(t  —  s,  z)E(s,  z)ds 


where  Q  is  determined  by  various  polarization  mechanisms  each  described  by  a  different  parameter  r, 
and  therefore  is  given  by 

g(t,z;F)  =  J  g(t,z\r)dF(r) 

where  T  C  [ti,T2].  In  particular,  if  the  distribution  were  discrete,  consisting  of  a  single  relaxation 
parameter,  then  we  would  again  have  (4). 

Note  that, 

P(t,  z)  —  /  Q(t  —  s,  z)E{s,  z)ds  +  £7(0,  z)E(t,  z)  +  £7(0,  z)E(t,  z)  (5) 

Jo 


where 

Q(t,  z;  F)  =  J  g(t,  z\  t )dF(r). 


Substituting  (5)  into  (1)  we  obtain 


fj,0eE(t,  z)  +  fj,0In  [<?  +  £7(0,  z)\  E(t,  z )  +  Mo-fa£7(0,  z)E(t,  z) 

+  [  Mo IciG(t— s,z)E(s,z)ds  — E"(t,z) 

Jo 

Further,  converting  to  weak  form  (and  multiplying  both  sides  by  c2),  we  find 


—Mo  Ja(t,  z). 


{erE(t,  •),  4>)  +  (-In  [a  +  67(0,  •)]  E(t,  •),  <f> >  +  <-/n£7(0,  -)E(t,  •),  J>) 
eo  eo 

+  (~  f  IaQ(t  -  s,  ■ )E(s ,  -)ds,  0)  -  ( c2E"{t ,  •),  J>)  =  -(-  J.(t,  ■),  </>), 
eo  Jo  eo 

where  er  =  ^  and  <j>  £  V  =  H] j(0,  1)  =  {</>  £  H1( 0,1)  :  0(1)  =  0}.  (Without  loss  of  generality,  we 
will  hereafter  assume  ir  =  1.)  Finally,  we  integrate  by  parts,  and  apply  the  boundary  conditions  (2)  to 
obtain 

(E(t,  •),  0)  +  {-In  [a  +  67(0,  2)]  E(t,  •),  0)  +  <- Jn£7(0,  -)E(t,  •),  0) 
eo  eo 

+  (-  f  InG(t  -  s,  ■ )E(s ,  -)ds,  0)  +  ( c2E\t ,  •),  0')  +  cE(t,  0)0(0)  =  •),  0), 

eo  Jo  eo 

which  we  can  rewrite  as 

(E,  0)  +  (7 E,  0)  +  </ 3E ,  0)  +  (  T  a(t  -  s,  ■ )E(s ,  -)ds,  0)  +  (c2E\  0')  +  cE(t,  0)0(0)  =  (J,<j>),  0  G  V,  (6) 

Jo 
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where 


7(2)  = 

=  —In  [a (z) 
eo 

+  G{  0,  z)] 

=  —In 
eo 

m- 

=  — InG(  0,  z 

)  =  -In 

[  9(0,  z 

eo 

eo 

It 

a(t,  z)  - 

=  —  InG(t,  z] 

)  =  -In  , 

[  g(t,z ; 

eo 

eo  J 

r 

—  J s  (£j  -2'. 

)• 

v(z)  +  J  g(0,z\T)dF(r) 


eo 


3  Estimation  Methodology 

Our  goal  is  to  estimate  the  probability  distribution  function  (PDF)  of  relaxation  parameters  F  £  V(T)  in 
a  given  model  of  the  polarization,  where  V(T)  is  the  set  of  all  PDFs  on  the  admissible  region  T  C  [n,  t2]. 
To  this  end  we  attempt  to  minimize  the  difference  between  model  simulations  and  observations  of  time- 
domain  data.  In  our  formulation  (for  details,  see  [BBL]),  the  observations,  Ej ,  are  of  the  electric  field 
E  at  discrete  times  tj  taken  at  z  =  0.  Each  simulation  is  a  solution  of  Maxwell’s  equation  given  in  (6) 
with  (3)  using  candidate  values  for  the  distribution  of  relaxation  parameters.  We  propose  a  standard 
least-squares  criterion  (which  is  equivalent  to  a  maximum  likelihood  estimation  in  many  situations-see 
[DG])  for  the  optimization  procedure  given  by 

J(F)  =  Y/\E(tj,0',F)-Ej\2  (7) 

3 

where  E(-,--F)  is  the  solution  of  (6)  with  (3)  corresponding  to  the  distribution  F.  Thus  the  inverse 
problem  is  to  solve 

min  J(F). 

FGV(T) 

In  practice,  one  may  choose  to  approximate  a  continuous  distribution  F  by  a  discrete  one  with,  for 
example,  N  elements.  This  approach  would  result  in  a  straight-forward  IV-dimensional  minimization 
problem.  Other  parameterizations  of  the  unknown  distribution  are  also  possible  to  reduce  the  problem 
to  a  finite  dimensional  one. 


4  Well-Posedness 

In  this  section  we  address  the  questions  of  well-posedness  of  the  form  of  Maxwell’s  equation  given  in  (6) 
with  respect  to  the  unknown  distribution  of  dielectric  parameters.  In  particular  we  wish  to  establish  the 
continuous  dependence  of  solutions  on  the  distributions  in  the  sense  of  the  Prohorov  metric.  This  leads 
in  turn  to  well-posedness  for  the  inverse  problems  involving  (7). 

First,  given  that  we  assume  g,  g,  and  g  are  uniformly  continuous  in  r  and  bounded  on  ([0,  T\  x  [0, 1]  x 
[n,T2])  where  0  <  n  <  r2  <  00,  then  Q,  Q,  and  Q  are  L°°  since,  for  example, 

\Q\00  <  J  |s(-,  ■;t)\00(IF(t)  <  Moo 

when  T  C  [ri,r2].  Therefore,  a ,  f3,  and  7  are  all  also  L°° .  Thus,  for  each  fixed  F  £  V(T),  a,  /3,  and  7 
satisfy  the  hypothesis  of  Theorem  1  in  [BBL,  p.  35],  so  that  given  $  £  V,  T  £  H  =  L2(0, 1),  a  unique 
solution  to  (6)  with  (3)  exists,  and  E  £  L2([0,T],V),  E  £  L2([0,  T],  H). 

It  remains  yet  to  show  the  continuous  dependence  of  solutions  on  F.  First  we  note  that  Fn  — ►  F  in 
the  Prohorov  sense  is  equivalent  to 

J  f{r)dFn(r)  — >  J  f{r)dF(r)  V  /GC[n,r2]. 

Since  g  is  uniformly  continuous  in  r,  we  have  then  that  a(t,z-,Fn)  — >  a(t,z;F)  a.e.  when  Fn  — >  F. 
We  will  define  a„  :=  a(t,Z',Fn)  and  a  :=  a(t,Z',F).  Because  |an| <  M  for  all  n,  the  Dominated 
Convergence  Theorem  implies  that  an  —>  a  in  L2([0,T]  x  [0, 1]).  Similarly,  we  can  show  that  (3n  — >  /3 
and  7„  — >  7  in  L2. 
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We  consider  arguments  for  fixed  (J  and  7.  We  need  to  show  that  (En,En)  — >  (E,  E)  in  L2([ 0,  T],V)x 
L2([0,  T\,  H)  when  an  —>  a,  where  En  denotes  the  solution  to  (6)  corresponding  to  a„. 

We  begin  by  subtracting  equation  (6)  corresponding  to  a  from  (6)  corresponding  to  a„  to  obtain 

(En  -  E,  c/>)  +  <7(iT  -  E),  <f>)  +  </ 3(En  -  E),  <f>)  +  (c2(En'  -  E'),  <f>') 

+  c  0)  —  E(t,  0)]  (f)(0)  =  —  (J  [an(t  —  s,  z)En(s,  z)  —  a(t  —  s,  z)E(s,  z)\  ds,  (f>). 

By  adding  and  subtracting  a  mixed  term  to  the  right  side,  we  have 

(En  —  E,<j>)  +  (7 (En  -  E),  4 >)  +  {P(En  -  E),  4>)  +  (c2(Enl  -  E'),  <f>') 

+  c[En(t,0)-E(t,0)^4>(0)  =  -{J^  an(En-E)ds,<f>)-{J^(an-a)Eds,4>).  (8) 

Following  the  general  procedure  in  [BBL]  (  see  also  [L],  [W]  for  related  fundamental  ideas  and  theory), 
we  approximate  En  and  E  in  Kn=span{u>i, ...,  wm}  by 

m 

Em(t,  z)  =  ^  e"(t)Wi(» 


Em(t,z )  =  y;  ej(t)wj(z). 


Then  (8)  becomes 


(EZ  -  Em,4>)  +  (7(S"  -  Em),  <f>)  +  (P(EZ  -  Em)  A)  +  (c(EZ  -  E’m)A) 

+  c  ^Em(t,  0)  -  Em(t,  0)j  4>(0)  =  —(j  a„(Em  —  Em)ds,4>)  ~  { J  («»  -  a)Emds,<f>), 

which  must  hold  for  all  <f>  £  Vm-  Since  Em  and  Em  are  both  in  Vm,  we  may  take  <f>  =  Em  —  Em.  Then 
we  have 

“  [l  Em-  Emil  +  cn  (EZ  -  Em,  £%  -  Em)]  +  |^7«  ~  Em) \l  +  c\Em(t,  0)  -  Emit,  0)|2 

=  (k(Em  —  Em),  Em  —  Em )  —  (  f  an(Em  —  Em)ds,  Em  —  Em)  —  (  f  («»  —  a)Emds,  Em  —  Em), 

Jo  Jo 

where  we  have  used  the  sesquilinear  form  <ti  :  V  x  V  — >  C  defined  by 

(Ji(4>A)  =  ( c24>'A')h  +  04>A)h 

for  </>,?/>  £  V  and  where  j3  =  (J  +  k  >  0  for  some  sufficiently  large  k  >  0. 

Integration,  along  with  the  Wellipticity  of  a  1  yields 

IfC(t)  -  Emit)  ll  +  Cl|  EZ(t)  -  Emit) \2V  +  2  IV7  K(«)  -  Em(s))\lds 

Jo 

+  2c\Em(t,  0)  -  Emit,  0)|2L2(Oit)  <  2  I  f*  FZiOdt I  ,  (9) 


FZiO  =  (k[E^(0  -  Emm,EUti)  ~  Emit))  -  (  -  Em(s)\ds,E^(0  -  Em(0) 

Jo 

-  (  f  [an  -  a\Em(s)ds,E £(£)  -  £m(0> 

Jo 

=  t2(0  +  t3(0  +  T4(0. 

Following  precisely  the  arguments  in  [BBL]  we  have  that 

fmm<  f  \\k2\Ei(o-Emmi  +  \\Ei(o-Emmi\  dt 
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and 


I T3(£M  <  Ki  /  | £%(()  -  Em(0\Ut  +  K2  /  | <(£)  -  Em(0\UC 


Lastly,  we  bound  the  final  term  by 


\nm  =  (  [an(Z-8)-a{Z-8)]Em(s)d8,E%l{Z)-Em{Z)) 


<o  /  {an(Z-s)-a(Z-s)]Em(S)ds  +  d<(0  -  £»(£)&• 


Since  Em  is  bounded  in  C([0,T],  V), 

,  i  /  r6  „  , 


\nm<n\  KB\an(H-s)-a(Z-s)\Hds)  +-|£”(0-Sm(0| 


<  /  |a„(£  -  a)  -  a(£  -  «)frdfl  +  -| £%(()  -  Em{t) \%. 


f  |r4(C)|dC  <  K3  f  ^  |a„(C  -a)-  a(£  -  a)\%dadi  +  Iu  f  \  i&(0  ~  Em($ \j 
Jo  Jo  Jo  Jo 


We  observe  that 


[  f  fit  -  s)dsd£  =  [  [  f(u)dsdu 
Jo  Jo  Jo  Jo 

=  /  f(u)(t  —  u)du 
Jo 


C  f  f(u)du  for  t,  u  €  [0,  T\. 
Jo 


Therefore 


|T4(CM  <  K3  /  \an(u)-a(u)\2Hdu  +  K4  |£&(0  -  Emit)\Ut- 


Using  the  above  bounds,  we  find  that  (9)  becomes 

\E^(t)-Em(t)\2H  +  Cl\E^(t)-Em(t)\2v  +  2  ^ \^(E^(s)-Em(s))\2Hds  +  2c\E^(;0)-Em(;0)\2L2i0tt) 

Jo 

<Ci  f  |0«)  -  Em(s)\2Hds  +  C2  [  |£”(s)  -  Em(s)\2Hds  +  C3  [  \an(s)  -  a(s)\2Hds.  (10) 
Jo  Jo  Jo 

Noting  that  |  •  \h  <  fJ>\  •  \v  for  some  constant  fj,  and  letting 

rT  2  2 

ACn  1=  /  |QJn(s)  —  Ol(s)|ifds  =  \otn  ~  Ot\ L2 ([0,T] ,H)  > 

Jo 

we  then  have  that  (10)  implies 

I Emit)  -  Emit) \2h  +  I £%(t)  -  Em(t)\2V  <  VKn  +  v  J*  [\E™  (s)  -  Em{s)\2v  +  | <(«)  -  Em{s) \2h]  ds 

for  some  v  >  0  independent  of  m  and  n.  Using  Gronwall’s  Inequality  we  have 

\E&{t)-Em(t)\%  +  \E&(t)-Em{t)fir<VKnevT  forfe[0,T].  (11) 

From  the  existence  proof  in  [BBL]  we  have  that  for  each  fixed  n,  Em  En  in  L2([0,  T],  V).  Therefore 
using  weak  lower  semicontinuity  of  norms  we  have 

\E  —  £'|i,2([o,t],v)  —  liminf  \Em  —  Em\L2([o,T] ,v)  > 
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with  a  similar  estimate  holding  for  En  —  E  in  the  L2{[ 0,  T\,  H)  norm.  We  may  apply  these  results  to  the 
integral  of  (11)  from  0  to  T  to  finally  obtain 

\En  —  £|i2([0,T],ir)  +  l-E1’1  —  -E|l2([o,t],v)  <  nKnTevT . 

Since  an  — i ►  a  in  L2,  then  nn  — >  0,  which  gives  continuous  dependence  of  ( E,E )  on  a.  Similar 
arguments  show  that  (E,  E)  depend  continuously  on  7  and  also  on  (3.  Thus  we  have  that  solutions  of 
(6)  with  (3)  depend  continuously  on  the  probability  measure  F  in  the  sense  that  the  map 

F  — >  (E,  E) 

is  continuous  from  V(T)  to  L2([0,T],V)  x  L2([0,  T],  H). 

Further,  this  yields  that  F  — >  J(F)  =  JT  \F(tj,0;  F)  —  Ej |2  is  continuous  from  V{T)  to  R1,  where 
V(T),  with  the  Prohorov  metric,  is  compact  for  T  compact.  Then  the  general  theory  of  Banks-Bihari  in 
[BB]  as  outlined  in  [BBPP]  can  be  employed  to  obtain  existence  and  stability  for  the  inverse  problem, 
as  well  as  an  approximation  theory  which  can  be  used  as  a  basis  for  a  computational  methodology. 


5  Conclusion 

We  have  presented  theoretical  results  on  a  model  for  the  electric  field  with  multiple  electric  polarization 
mechanisms  in  a  dielectric  material.  This  provides  a  firm  foundation  for  an  inverse  problem  formulation 
to  determine  an  unknown  probability  distribution  of  parameters  which  describe  the  dielectric  properties 
of  the  material.  To  this  end,  we  have  shown  the  continuous  dependence  of  the  solutions  with  respect 
to  the  unknown  distributions  in  the  Prohorov  metic.  This  argument,  combined  with  previous  results  on 
existence  and  uniqueness  in  Maxwell  systems,  demonstrate  the  well-posedness  of  the  model.  Moreover, 
the  theory  described  in  [BB] ,  [BBPP]  can  be  combined  with  our  results  here  to  provide  existence,  stability, 
and  an  approximation  theory  for  the  associated  inverse  problems. 


Acknowledgements 

This  research  was  supported  in  part  by  the  U.S.  Air  Force  Office  of  Scientific  Research  under  grant 
AFOSR  F49620-01-1-0026  and  in  part  by  the  NASA  Langley  GSRP  under  grant  NGT-1-02004. 


References 

[BB]  H.  T.  Banks  and  K.  L.  Bihari,  Modeling  and  estimating  uncertainty  in  parameter  estimation, 
Inverse  Problems,  17(2001),  pp.  1-17. 

[BBPP]  H.  T.  Banks,  D.  M.  Bortz,  G.  A.  Pinter,  L.  K.  Potter,  Modeling  and  Imaging  Techniques 
with  Potential  for  Application  in  Bioterrorism,  Tech.  Rep.  CRSC-TR03-02,  Center  for  Re¬ 
search  in  Scientific  Computation,  North  Carolina  State  University,  January,  2003;  Chapter 
6  in  Bioterrorism:  Mathematical  Modeling  Applications  in  Homeland  Security,  (H.T.  Banks 
and  C.  Castillo-Chavez,  eds.),  SIAM  Frontiers  in  Applied  Math,  FR28,  pp.  129-154,  SIAM, 
Philadelphia,  2003. 

[BBL]  H.  T.  Banks,  M.  W.  Buksas,  and  T.  Lin,  Electromagnetic  Material  Interrogation  Using  Conduc¬ 
tive  Interfaces  and  Acoustic  Wavefronts,  SIAM  Frontiers  in  Applied  Mathematics,  Philadelphia, 
2000. 

[B]  P.  Billingsley,  Convergence  of  Probability  Measures,  Wiley,  New  York,  1968. 

[DC]  M.  Davidian  and  D.  Giltinan,  Nonlinear  Models  for  Repeated  Measurement  Data,  Chapman 
and  Hall,  London,  1998. 

[E]  R.  S.  Elliott,  Electromagnetics:  History,  Theory  and  Applications,  IEEE  Press,  New  York, 

1993. 

[L]  .1.  L.  Lions,  Optimal  Control  of  Systems  Governed  by  Partial  Differential  Equations,  Springer- 

Verlag,  New  York,  1971. 

[W]  .1.  Wloka,  Partial  Differential  Equations,  Cambridge  University  Press,  1992. 


6 


